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Abstract 

Today, the human brain can be studied as 
a whole. Electroencephalography, magnetoen- 
cephalography, or functional magnetic resonance 
imaging (fMRI) techniques provide functional 
connectivity patterns between different brain ar- 
eas, and during different pathological and cog- 
nitive neuro-dynamical states. In this Tutorial 
we review novel complex networks approaches to 
unveil how brain networks can efficiently man- 
age local processing and global integration for 
the transfer of information, while being at the 



same time capable of adapting to satisfy chang- 
ing neural demands. 

1 Introduction 

In recent years, complex networks have provided 
an increasingly challenging framework for the 
study of collective behaviors in complex systems, 
based on the interplay between the wiring ar- 
chitecture and the dynamical properties of the 
coupled units [U [2]. Many real networks were 
found to exhibit small- world features. Small- 
world (SW) networks are characterized by hav- 



ing a small average distance between any two 
nodes, as random graphs, and a high clustering 
coefficient, as regular lattices [3J, SI El [6] . Thus, a 
SW architecture is an attractive model for brain 
connectivity because it leads distributed neural 
assemblies to be integrated into a coherent pro- 
cess with an optimized wiring cost [H H2 [9] . 

Another property observed in many networks 
is the existence of a modular organization in the 
wiring structure. Examples range from RNA 
structures, to biological organisms and social 
groups. A module is currently defined subset 
of units within a network such that connections 
between them are denser than connections with 
the rest of the network. It is generally acknowl- 
edged that modularity increases robustness, flex- 
ibility and stability of biological systems [TOl [11] . 
The widespread character of modular architec- 
ture in real-world networks suggests that a net- 
work's function is strongly ruled by the organi- 
zation of their structural subgroups. 

Recent studies have attempted to character- 
ize the functional connectivity (patterns of sta- 
tistical dependencies) observed between brain 
activities recorded by electroencephalography 
(EEG), magnetoencephalography (MEG), or 
functional magnetic resonance imaging (fMRI) 
techniques [121 El El El E]- Surprisingly, 
functional connectivity patterns obtained from 
MEG and EEG signals during different patholog- 
ical and cognitive neuro-dynamical states, were 
found to display S W attributes [El E| ; whereas 
functional patterns of fMRI often display a struc- 
ture formed by highly connected hubs, yielding 
an exponentially truncated power law in the de- 
gree distribution [121 El El- For a complete 
review of these issues, reader can refer to the 
Refs. PHIig. 

In functional networks, two different nodes 
(representing two electrodes, voxels or source 



regions) are supposed to be linked if some de- 
fined statistical relation exceeds a threshold. 
Regardless of the modality of recording activ- 
ity (EEG, MEG or fMRI), topological features 
of functional brain networks are currently de- 
fined over long periods of time, neglecting pos- 
sible instantaneous time-varying properties of 
the topologies. Nevertheless, evidence suggests 
that the emergence of a unified neural process is 
mediated by the continuous formation and de- 
struction of functional links over multiple time 
scales [201 El EI]. 

Empirical studies have lead to the hypothesis 
that transient synchronization between distant 
and specific neural populations underlies the in- 
tegration of neural activities as unified and co- 
herent brain functions [19]. Specialized brain 
regions would be largely distributed and linked 
to form a dynamical web-like structure of the 
brain [20] . Thus, brain regions would be par- 
titioned into a collection of modules, represent- 
ing functional units, separable from -but related 
to- other modules. Characterizing the dynami- 
cal modular structure of the brain may be crucial 
to understand its organization during different 
pathological or cognitive states. An important 
question is whether the modular structure has 
a functional role on brain processes such as the 
ongoing awareness of sensory stimuli or percep- 
tion. 

To find the brain areas involved in a given 
cognitive task, clustering is a classical approach 
that takes into account the properties of the 
neurophysiological time series. Previous stud- 
ies over the mammalian and human brain net- 
works have successfully used different methods 
to identify clusters of brain activities. Some clas- 
sical approaches, such as those based on princi- 
pal components analysis (PCA) and independent 
components analysis (ICA), make very strong 
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statistical assumptions (orthogonality and sta- 
tistical independence of the retrieved compo- 
nents, respectively) with no physiological justifi- 
cation [221 125) . 

In this Tutorial, we review an approach that 
allows to characterize the dynamic evolution of 
functional brain networks \24\ [25] . We illustrate 
this approach on connectivity patterns extracted 
from MEG data recorded during a visual stimu- 
lus paradigm. Results reveal that the brain con- 
nectivity patterns vary with time and frequency, 
while maintaining a small-world structure. Fur- 
ther, we are able to reveal a non-random mod- 
ular organization of brain networks with a func- 
tional significance of the retrieved modules. This 
modular configuration might play a key role in 
the integration of large scale brain activity, fa- 
cilitating the coordination of specialized brain 
systems during a cognitive brain process. 



2 Materials and methods 

To illustrate our approach, we consider the brain 
responses recorded during the visual presenta- 
tion of non-familiar pictures. Although our ap- 
proach is applicable to any of the functional 
methods available (EEG, fMRI, MEG), here we 
use the magnetoencephalography. This modal- 
ity of acquisition has the major feature that col- 
lective neural behaviors, as synchronization of 
large and sparsely distributed cortical assem- 
blies, are reflected as interactions between MEG 
signals [26]. We study the functional connec- 
tivity patterns associated with dynamic brain 
processes elicited by the repetitive application 
(trials) of a external visual stimulus [27]. For 
this experiment, a collection of 48 simple struc- 
tural images and scrambled images were ran- 
domly shown to epileptic patients for a peri- 



ode of 150 ms with an inter-stimulus interval 
of 2 s. Patients were required to respond by 
pressing a button each time an image was per- 
ceived. The event-related brain responses were 
recorded (from two patients) with a whole-head 
MEG system (151 sensors; VSM MedTech, Co- 
quitlam, BC, Canada) digitized at 1.25 kHz with 
a bandpass of — 200 Hz. 

The basic steps of our approach are schemat- 
ically illustrated in Fig. [TJ Each of the sig- 
nals is decomposed into time-frequency compo- 
nents, as shown in panel a). The relations be- 
tween two signals j and k are firstly defined 
in time- frequency space, as shown in panel b). 
A statistical criterion is then used to define a 
functional connectivity matrix for each time- 
frequency point, panel c). The details of the 
statistical criterion we adopted are reported in 



Section 2.1 In panel d), topological metrics are 
extracted from the connectivity patterns to ob- 
tain a time-frequency characterization of brain 
networks. The metrics investigated are analyzed 
in Section 2.2 Finally, in panel e), at a given 
frequency, or time instant of interest, the modu- 
lar structure is characterized as discussed in Sec- 
tions 12.31 and 12.41 To evaluate the features of 
brain connectivity, the obtained functional net- 
works are compared with equivalent regular and 
random networks. 

2.1 Estimation of functional connec- 
tivity 

A unified definition of brain connectivity is diffi- 
cult from the fact that the recorded dynamics re- 
flect the activities of neural networks at different 
spatial and temporal resolutions. Three types of 
connectivity are currently considered: anatomi- 
cal (description of the physical connections be- 
tween two brain sites), functional (defined by a 
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Figure 1: General scheme for the extraction of the time- varying brain networks: (a) signals are 
decomposed into time-frequency components to compute (b) pair-wise relations; (c) functional 
connectivity matrices are extracted at each point of the time-frequency space, defining (d) the 
functional brain networks used to extract the topological attributes (color codes the nodes degree) 
and the (e) modular structure (brain sites belonging to each module are arbitrarily colored). See 
details in the text. 
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temporal correlation between distant neurophys- 
iological events) and effective (causal influence 
that a neural system may exert over another). 
Here, we consider the functional links in brain 
signals defined by means of the phase-locking 
value (PLV) computed between all pairs of sen- 
sors [28]. To compute the PLV values, we used 
a complex Morlet's wavelet function defined as 
w(t,f ) = Aexp(-t 2 /2a?) x ex.p(i2irf t). Nor- 
malization factor A was set to A = (cxiy^) -1 / 2 . 
at = m/27r/o, m is a constant that defines the 
compromise between time and frequency reso- 
lution, and /o is the center frequency of the 
wavelet. Hence, in time domain, its real and 
imaginary parts are a cosine and a sine, re- 
spectively, of which the amplitude envelope is 
a Gaussian with a standard deviation of at- In 
frequency domain, the Morlet wavelet is also a 
Gaussian with a standard deviation cry given 
m = fo/crf. Here, m was chosen to be 7. By 
means of this complex wavelet transform an in- 
stantaneous phase (j>l rial (t, f) is obtained for each 
frequency component of signals i = 1, . . . , M at 
each repetition of the stimulus (trial). The PLV 
between any pair of signals (i, k) is inversely re- 
lated to the variability of phase differences across 
trials: 



PLV i)fc (t,/) 



1 



V/ 



trials 



N f , 



exp J 



l (t,f)-4>p al (t,f)) 



trial=l 

where N tr i a is is the total number of trials. If 
the phase difference varies little across trials, its 
distribution is concentrated around a preferred 
value and PLV ~ 1. In contrast, under the null 
hypothesis of a uniformity of phase distribution, 
PLV values are close to zero. 

Finally, to assess whether two different sen- 
sors are functionally connected, we calculated 



the significance probability of the PLV values 
by a Rayleigh test of uniformity of phase. Ac- 
cording to this test, the significance of a PLV 
value determined from N tr iais can be calculated 
as p = exp(— NtriaisP^N 2 ) [29] • To correct for 
multiple testing, the False Discovery Rate (FDR) 
method was applied to each matrix of PLV val- 
ues |30j. With this approach, the threshold of 
significance PLVth was set such that the ex- 
pected fraction of false positives is restricted to 
q < 0.05. 

In the construction of the networks, a func- 
tional connection between two brain sites was 
assumed as an undirected and unweighted edge 
(Aij = 1 if PLVij > PLVth] and zero other- 
wise). Although topological features can also 
be straightforwardly generalized to weighted net- 
works [32], we obtained qualitative similar re- 
sults (not reported here) for weighted networks 
with a functional connectivity strength between 
nodes given by Wij = PLVij. More refined sta- 
tistical tools can also be used to estimate time- 
varying and directed brain networks [31] . 

2.2 Time- varying structure of brain 
networks 

A set of metrics can be used to characterize the 
topological properties of the functional networks 
we have constructed [U [2]. Here, we use three 
key parameters: mean degree (K), clustering in- 
dex C and global efficiency E. Briefly, the degree 
ki of node i denotes the number of functional 
links incident with the node and the mean de- 
gree is obtained by averaging ki across all nodes 
of the network. The clustering index quantifies 
the local density of connections in a node's neigh- 
borhood. The clustering coefficient Cj of a node 
i is calculated as the number of links between 
the node's neighbors divided by all their pos- 
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sible connections and C is defined as the aver- 
age of Cj taken over all nodes of the network [3J. 
The global efficiency E provides a measure of 
the network's capability for information transfer 
between nodes and is defined as the inverse of 
the harmonic mean of the shortest path length 
Lij between each pair of nodes [6]. The node- 
efficiency Ei of the i th node is likewise defined 
as the inverse of the harmonic mean of the min- 
imum path length between node i and all other 
nodes in the network. 

To asses the small- world behavior of functional 
networks, we perform a benchmark comparison 
of the functional connectivity patterns [3J. For 
this, the clustering and efficiency coefficients of 
functional networks are compared with those ob- 
tained from equivalent random and regular con- 
figurations. Regular networks were obtained by 
rewiring the links of each node to its nearest (in 
the sensors space) neighbors, yielding a nearest- 
neighbor connectivity with the same degree dis- 
tribution as the original network. To create an 
ensemble of equivalent random networks we use 
the algorithm described in Ref . [3j . According to 
this procedure, each edge of the original network 
is randomly rewired avoiding self and duplicate 
connections. The obtained randomized networks 
preserve thus the same mean degree as the origi- 
nal network whereas the rest of the wiring struc- 
ture is random. 

2.3 Network modularity 

Many real networks have a modular structure, 
i.e. their associated graphs are in general glob- 
ally sparse but locally dense. In these net- 
works, modules are defined as groups of vertices 
linked such that connections between them are 
denser than connections with the rest of the net- 
work. It is currently accepted that a partition 



V = {C\, . . . ,Cm} represents a good division in 
modules if the portion of edges inside each mod- 
ule C{ (intra-modular edges) is high compared 
to the portion of edges between them (inter- 
modular edges). The modularity Q(V), for a 
given partition V of a network is formally de- 
fined as |33j: 



M 



Q{V) = 



Is ( ks \ 
L ~ \2LJ 



(1) 



where M is the number of modules, L is the 
total number of connections in the network, l s 
is the number of connections between vertices in 
module s, and k s is the sum of the degrees of the 
vertices in module s. 

To partition the functional networks in mod- 
ules, we used a random walk-based algo- 
rithm [33], because of its ability to manage 
very large networks, and its good performances 
in benchmark tests [34|, 13 7j . Similar theoreti- 
cal frameworks have been recently proposed for 
spectral coarse-graining [35| [36] . The algorithm 
is based on the intuition that a random walker on 
a graph tends to remain into densely connected 
subsets corresponding to modules. Let Pij = 
to be the transition probability from node i to 
node j, where Aij denotes the adjacency matrix 
and k{ is the degree of the i th node. This defines 
the transition matrix {P l )ij for a random walk 
process of length t (denoted here P/- for simplic- 
ity). One can notice that, if two vertices i and j 
are in the same community, the probability P\- 



is high, and P\ k 

The metric used to quantify the structural 
similarity between vertices is given by 



Pij 



N 

l=i 



(Pfl 



pt \2 



(2) 
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This distance has several advantages: it quanti- 
fies the structural similarity between vertices and 
it can be used in an efficient clustering algorithm 
to maximize the network modularity Q. Fur- 
ther, using matrix identities, the distance p can 
be written as pfj = YZ=2 ( v a(i) ~ v a (j)) 2 ; 
where (\ a )i^a^n and (v a )i^ a ^n are the n eigen- 
values and right eigenvectors of the matrix P, re- 
spectively [34 . This relates the random walk al- 
gorithm to current methods using spectral prop- 
erties of the graphs [6S\ I36| . The random- walk 
based approach, however, needs not to explicitly 
compute the eigenvectors of the matrix; a com- 
putation that rapidly becomes intractable when 
the size of the graphs exceeds some thousands of 
vertices. 

To find the modular structure, the algorithm 
starts with a partition in which each node in 
the network is the sole member of a module. 
Modules are then merged by an agglomera- 
tive approach based on a hierarchical clustering 
method [39]. The algorithm stops when all the 
nodes are grouped into a single component. At 
each step the algorithm evaluates the quality of 
partition Q. The partition that maximizes Q is 
considered as the partition that better captures 
the modular structure of the network. In the cal- 
culation of Q, the algorithm excludes small iso- 
lated groups of connected vertices without any 
links to the main network. However, these iso- 
lated modules are considered here as part of the 
network for the calculation of the topological pa- 
rameters. 

2.4 Comparison of network partitions 

To evaluate the agreement between modules as- 
signments at a given time instant or frequency 
one can used the adjusted Rand index Ra [IT] , 
which is a traditional criterion for comparison of 



different results provided by classifiers and clus- 
tering algorithms, including partitions with dif- 
ferent numbers of classes or clusters. For two 
partitions P and P' , the original Rand index is 
defined as [10] R = a+ %^ +d ; where a is num- 
ber of pairs of data objects belonging to the 
same class in P and to the same class in P' , 
b is number of pairs of data objects belonging to 
the same class in P and to different classes in 
P' , c is the number of pairs of data objects be- 
longing to different classes in P and to the same 
class in P', and d is number of pairs of data ob- 
jects belonging to different classes in P and to 
different classes in P' . Thus the index R has 
a straightforward interpretation as a percentage 
of agreement between the two partitions and it 
yields values between (if the two partitions are 
randomly drawn) and 1 (for identical partition 
structures) . 

The Rand index, however, has a bias if a par- 
tition is composed by many clusters, and it can 
take a non-null value for two completely ran- 
dom partitions. The index R can be straight- 
forwardly corrected for the expected value un- 
der the null hypothesis according to the follow- 
ing general scheme: Ra = — '^L,. Us- 

° ° max{ -ft} — E{H\ 

ing the generalized hypergeometric distribution 
as the null hypothesis, the adjusted Rand index 
that corrects for the expected number of nodes 
pairs placed in the same module under two ran- 
dom partitions is given by |41j 

_ (a+c)(q+6) 

Ra = a+b+c+d (3) 

2a+b+c _ (a+c)(a+b) v 1 
a+b+c+d a+b+c+d 

which has an expected value of zero under the 
null hypothesis, and it takes a maximum value of 
one for a perfect agreement of the two partitions. 
Thus, the adjusted Rand index is a statistics on 
the level of agreement or correlation between two 
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partitions. 

3 Results 

3.1 Time-frequency dependence of 
brain networks 

Fig. [2] shows the topological attributes of func- 
tional networks elicited by the -unexpected- im- 
ages. Pictures show the values of the mean de- 
gree, clustering index and efficiency of networks 
between, calculated at each point of the time- 
frequency space, 600 ms before and 1 s after the 
onset of the stimulus. 

The first crucial observation is that functional 
connectivity patterns are not time- invariant, but 
instead they exhibit a rich time-frequency struc- 
ture during the neural processing. All the 
topological features (specially (K) and C) ex- 
hibit high values in a frequency band close to 
10 Hz, which is a spectral component mostly 
involved in the processing of visual informa- 
tion [27] . Whereas the functional networks in the 
frequency range of 10 — 30 Hz display large pat- 
terns of synchronization/desynchronization be- 
fore the stimuli, a highly connected pattern is 
induced by the stimulus at about 250 ms and 
between 15 and 25 Hz, suggesting a connectiv- 
ity induced by the unexpected sensory stimuli. 
This is followed by weak connected structures 
at frequency bands close to 7 and 15 Hz arising 
during the post-stimulus activities and marking 
the transition between the moment of percep- 
tion and the motor response of the subject. The 
topological features of these connectivity pat- 
terns were detected as statistically different from 
the pre-stimulus epoch by a Z-test corrected by 
a FDR at q < 0.05. Brain activities above 30 Hz 
are characterized by a poor global connectivity. 
Local parameters, ki, c% and Ei, for each sensor 



of the network are shown at three different time 
instants for a frequency of 20 Hz. During the 
processing of the stimulus, a time-space variabil- 
ity of connectivity is observed. Before the onset 
of the stimulus, the networks are characterized 
by a very sparse connectivity. Then, a clear clus- 
tered structure triggered by the stimulus appears 
at t = 250 ms, defining two main regions (frontal 
and occipital) with a high density of connections. 
After the stimulus, the functional wiring displays 
again a sparse structure. 

3.2 Small-world behavior of brain net- 
works 

The comparison of the brain networks against 
random and regular configurations is shown 
Fig. [3j Typically, small- world networks exhibit a 
E sw greater than regular lattices, but less than 
random wirings E[ at < E sw < E rnc i; while for 
the mean cluster index, C rn d < C sw < C\ a t is ex- 
pected [3]. Results reveal that, despite the vari- 
ability observed, functional networks display a 
topology different from regular and random net- 
works. Namely, /r F > > 1 and ,J? > < 1, which 
indicates a SW structure ((...) stays for an aver- 
age over the ensembles of equivalent networks). 

(E ) (E ) 

Further, v l £ u < 1 and v r £ d/ > 1, supporting 

the hypothesis of a SW connectivity. 

It is important to emphasize that, in con- 
trast with previous studies which have focused 
on time-invariant networks |12[ [TTfl [T^l [T5"l 116] . 
our approach reveals a dynamical small-world 
connectivity at multiple time scales. This is a 
remarkable result, insofar as it suggests that the 
processing of a stimulus involves an optimized 
(in a SW sense) functional integration of distant 
brain regions by a dynamic reconfiguration of 
links. 
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Figure 2: Time-frequency maps of topological features extracted from brain networks associated 
to a visual stimulus presentation (arriving at t = 0). (a) mean degree {K), (b) clustering index C 
and (c) efficiency E. The reported values refer to the average over subjects. Dotted lines outline 
the regions revealing a significant change from the pre-stimulus region. Lower row: topographic 
distribution of the local parameters for the 20 Hz activities (indicated by the thick dashed line) at 
three different time instants. 

ijgi: mm 

£*M if: ^^Hi 

■0.5 0.5 1 -0-5 0.5 1 

time (s) time (s) 

Figure 3: Comparison of functional networks with random and regular configurations: time- 
frequency maps of (a) C/{Ci a t), (b) C/(C rn d), (c) (E[ at }/E and (d) (E rn d) /E. Results of equivalent 
random and regular networks refer to the average of 20 realizations. 

3.3 Evolution of functional modules 




0.5 

time (s) 



time (s) 



A potential modularity of brain-webs is sug- 
gested by the fact that brain networks display a 



clustering index larger than that obtained from 
random configurations [42J. Indeed, the presence 
of modules is actually confirmed by the high val- 
ues of Q obtained for brain networks extracted 
from brain activities at different time instants 
and frequencies. Fig. [4] shows the spatial dis- 
tribution of the modules for different networks 
with the following Q values: (a) Q = 0.55, (b) 
Q = 0.33, (c) Q = 0.53, (d) Q = 0.53, (f) 
Q = 0.49 and (g) Q = 0.53. 

Results show that brain networks have a time- 
varying structure with a number of modules that 
changes with time and frequency. From plots 
one can observe that, despite the spatial variabil- 
ity observed at different time instants and fre- 
quencies, functional modules fit well some known 
brain regions including visual, somatosensory 
and auditory processing areas. Before the on- 
set of the stimulus, however, networks are char- 
acterized by modular structures that define four 
main regions (anterior and posterior for both left 
and right hemisphere). We notice that the par- 
titions of both networks present a relatively low 
agreement yielding a Ra = 0.39. Then, the large 
connectivity triggered by the stimulus is also ac- 
companied by an increase in the number of mod- 
ules, yielding a more complex modular struc- 
ture. The observed changes are directly related 
with the specific nature of the task:the detection 
and low-level processing of the stimulus involves 
the visual system, but further processing as the 
identification and perception of the picture re- 
quires the mediation of regions as those located 
in frontal regions. Surprisingly antero-posterior 
relations elicit a large and unique module fitting 
fronto-occipital regions. Although a one-to-one 
assignment of anatomo-functional roles to each 
detected module is difficult to define, results re- 
veal some other interesting modules, as the ones 
located over the motor cortex at / = 18.5 Hz. 



Then the post-stimulus activities recovers again 
a simpler spatial organization of modules. It is 
worthy to notice that the pre- and post-stimulus 
networks have a very similar modular architec- 
ture only for the brain activities at the frequency 
band of 10 Hz. This high agreement is con- 
firmed by a high value of the adjusted Rand in- 
dex (Ra = 0.626), compared with the values less 
than 0.38 obtained for other frequencies. 

These are remarkable results as they support 
the hypothesis that brain dynamics relies on dif- 
ferent modular organizations to integrate distant 
specialized, but functionally related, brain re- 
gions. Our findings suggest modularity as an 
organization basis leading distributed groups of 
specialized neural assemblies to be integrated 
into a coherent process during different cognitive 
or pathological states. A modular description of 
brain networks might provide, more in general, 
meaningful insights into the functional organiza- 
tion of brain activities during others neural func- 
tions, such as attention and consciousness. 

4 Conclusion 

In this Tutorial we have addressed a fundamen- 
tal problem in brain networks research: whether 
and how brain behavior relies on the coordination 
of a dynamic mosaic of functional brain mod- 
ules during cognitive states. We have proposed a 
method to study the time-frequency dependen- 
cies of functional brain networks, thus offering 
an instantaneous description of the brain archi- 
tecture. Applied to a visual stimulus paradigm, 
the method reveals that the functional brain con- 
nectivity evolves in a small-world structure dur- 
ing the different episodes of the neural process- 
ing. Furthermore, by using a random walk-based 
analysis, we have identified a non-random modu- 
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(e) (f) (g) 



Figure 4: Topographical distribution of the modules extracted from brain networks at different 
time instants and frequencies: (a) time instant t = -0.25 s, frequency f =20 Hz; (b) t = 0.25 s after 
the presentation of the stimulus at f =20 Hz; (c) t =0.75 s, f =20 Hz; (e) t =-0.25 s, f =10 Hz; (f) 
t =0.25 s, f =10 Hz and (g) t =0.75 s, f =10 Hz. Brain sites belonging to each functional brain 
module were arbitrarily colored (there is no color correspondence between the modules of different 
networks). For the sake of clarity, isolated nodes were colored in black, (d) Time- frequency maps 
of mean degree is plotted to help network's localization in the time frequency space 
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lar structure in the functional brain connectivity. 

The present analysis was performed on MEG 
data in sensor space, which contains some in- 
herent spurious correlation between magnetic 
fields on the surface of the brain. Although 
this caveat does not affect the characterization 
of the global network topology, accurate infer- 
ences about anatomical locations needs a source 
reconstruction of the activity in the cortex. In 
this study, we have reduced the influence of spu- 
rious correlations by simply excluding the near- 
est sensors from the computation of PLV values. 

Our approach may provide meaningful in- 
sights into how brain networks can efficiently 
manage a local processing and a global integra- 
tion for the transfer of information, while being 
at the same time capable of adapting to satisfy 
changing neural demands. Although the neuro- 
physiological mechanisms involved in the func- 
tional integration of distant brain regions are 
still largely unknown, a dynamic SW organiza- 
tion is a plausible solution to the apparently op- 
posing needs of local specificity of activity ver- 
sus the constraints imposed by the coordination 
of distributed brain areas. The modular struc- 
ture constitutes therefore an attractive model for 
the brain organization as it supports the coexis- 
tence of a functional segregation of distant spe- 
cialized areas and their integration during brain 
states [H [5] . We suggest that this network de- 
scription might provide new insights into the un- 
derstanding of human brain connectivity during 
pathological or cognitive states. 

Applied to other multivariate data, our ap- 
proach could provide new insights into the struc- 
ture of the time- varying connectivity at a certain 
time [24] . A modular description of brain net- 
works might provide, more in general, meaning- 
ful insights into the functional organization of 
brain activities recorded with others neuroimag- 



ing techniques (EEG, MEG or fMRI) during di- 
verse cognitive or pathological states [23]. In 
this study, the functional links have been defined 
in MEG signals by means of the phase-locking 
value. We notice, however, that other time- 
frequency methods (e.g. wavelet cross-spectra) 
can also be used to detect and characterize a 
time-varying connectivity of spatially extended, 
nonstationary systems (e.g. financial or epidemi- 
ological networks). 
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